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Modified  Hertz  theory  for  spherical  indentation 


E.H.  VToff.' 

Department  of  Metallurgy  and  Materials  Science, 
University  of  Cambridge,  Cambridge  CB2  3QZ,  England 


Abstract 

The  elastic  field  in  a  half-space  indented  by  a  rigid  sphere  is 
accurately  described  by  Hertz1  equations  if  the  contact  radius  is  small 
compared  with  that  of  the  sphere.  For  wider  contact  the  equations  are 
not  expected  to  be  reliable,  but  are  nevertheless  often  used  for 
convenience.  This  paper  describes  some  of  the  errors  which  may  develop 
as  contact  widens,  and  presents  a  second  elastic  field  which  provides 
a  first  order  correction.  From  this  it  is  possible  to  estimate  the 
probable  error  introduced  in  a  particular  variable  by  applying  the  Hertz 
equations  beyond  the  recommended  range . 

Some  useful  relations  between  sane  familiar  elastic  indentation 


fields  are  also  reported. 


1. 


Introduction 


It  was  clearly  stated  by  Hertz  (1881)  that  his  solution  for  the 
elastic  field  due  to  pressure  between  spheres  should  be  applied  only  to 
contact  regions  of  radius  small  compared  with  that  of  the  spheres.  His 
formulae  have,  for  many  years,  been  applied  to  the  indenting  of  a  flat 
elastic  specimen  by  a  hard  sphere,  and  found  accurate  and  reliable  in  the 
range  a/R  <  0.1,  where  a,  R  are  the  radii  of  contact  and  of  the  sphere 
respectively. 

The  theory  has  also  been  used  beyond  this  value,  and  the  Hertz 
formulae  now  form  the  basis  of  an  ever-widening  range  of  methods  for 
testing  and  determining  the  properties  of  materials.  Because  of  this,  it 
is  of  interest  to  calculate  another  term  in  the  elastic  solution,  extending 
its  range  to  larger  a/R,  in  order  to  make  clear  the  conditions  under  which 
this  second  term  may  safely  be  ignored. 

In  this  paper  we  imagine  a  rigid  sphere  pressed  on  an  elastic 
half-space,  and  calculate  the  displaced  form  of  the  initially  flat  surface 
using  Hertz'  equations.  It  is  shown  that  there  is  some  mi'Sfit  between  the 
deformed  surface  and  the  sphere,  and  a  second  elastic  field  is  developed  as 
a  correction.  The  theory  of  this  second  stress  field  is  described  in 
detail  in  $3,  and  some  relevant  comments  of  a  more  general  nature  are  made 
in  S4.  The  latter  section  may,  however,  be  omitted  by  the  reader  who  wishes 
to  proceed  to  the  applications  in  §5.  Finally  the  method  of  correction  is 
demonstrated  in  the  practical  case  of  determining  the  elastic  modulus  E  of 
a  pliant  material. 

The  displacements,  stresses  and  strains  are  those  of  first  order 
linear  elasticity  theory  throughout. 


For  a  rigid  spherical  indenter  pressed  with  load  F  on  the  flat 


surface  of  an  elastic  half-space, 
as  given  by  the  Hertz  theory  are: 


the  components  of  displacement  u  and  w 
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within  the  contact  area  o  <  a.  Here  v  and  G  are  Poisson's  ratio  and  the 
shear  modulus  respectively,  and  the  cylindrical  axes  p  are  as 

indicated  in  fig.l. 

Near  the  2  axis,  these  displacements  give  very  close  fit  of  the 
deformed  surface  to  a  sphere  of  radius  R  ,  where 

1/Rh  =  3(1  -  v) F/8G  aJ  (2) 

since  the  w  displacement  defines  a  parabola  with  this  value  of  maximum 
curvature.  But  as  p  increases  the  parabola  diverges  from  the  sphere,  and 
the  w  displacement  therefore  introduces  a  growing  gap,  while  the  u 

P 

component  tends  to  close  it.  For  large  values  of  Poisson's  ratio, v,  the 
component  w  has  the  greater  effect,  since  up  -*■  0  as  v  -*•  0.5.  So  for  large 
v  the  Hertz  solution  creates  a  gap  between  sphere  and  surface  in  the  outer 
part  of  the  contact  region,  despite  the  transmission  of  pressure  there. 

In  other  words,  contact  is  lost  over  part  of  the  "contact  area”. 

For  small  values  of  v  the  component  up  is  dominant,  bringing  the 
indented  surface  actually  within  the  volume  of  the  sphere.  For  intermediate 
values  of  v  there  is  an  approximate  balance  between  the  two,  gap  or 
overlap  depending  on  the  extent  of  contact  a/R. 

In  fig. 2  are  curves  indicating  the  gaps  and  overlaps  implicit  in  the 


calculated  at  p  =  0.9a  in  each  case,  since  this  is  near  the  position  of 

maximum  u  .  The  error  for  small  v  is  always  an  overlap,  which  increases 

0 

with  decreasing  v  and  with  increasing  a/R. 

These  errors,  although  small,  show  a  logical  inconsistency  in 
applying  the  Hertz  theory  outside  the  range  a/R  .2,  since  the  prescribed 
pressure  cannot  readily  be  transmitted  across  a  gap,  nor  can  sphere  and 
specimen  occupy  the  same  point  in  space.  The  value  of  the  contact  radius, 
a,  for  a  given  load  and  sphere  radius,  becomes  doubtful  also,  but  it  is 
not  immediately  apparent  whether  larger  or  smaller  values  would  reduce  the 
errors  in  fig. 2.  Since  the  values  of  a  are  commonly  used  in  estimates  of 
strength  and  toughness  derived  from  the  indentation  of  hard,  brittle 
materials,  the  question  of  their  accuracy  becomes  important. 

To  investigate  this,  better  fitting  contact  may  be  made  by  combining 
the  Hertz  field  with  another  exact  elastic  one  (the  parabolic  pressure- 
field)  in  proportions  1  +  dt*d,  where  d  is  a  small  parameter,  so  that 
each  of  the  two  separate  solutions,  and  the  combined  one,  applies  the  same 
load  F  to  the  same  contact  area  p  <  a.  This  means  that  a  is  fixed,  while 
the  parameter  d  is  varied  and  adjusted  so  as  to  make  the  displaced  surface 
as  nearly  spherical  as  possible  over  the  whole  contact  area.  The  value  of 
d  is  found  in  this  way  for  each  value  of  F,  and  thence  a  modified  form  of 
equation  (2) . 

But  first  the  new  elastic  field  must  be  described. 

3.  The  parabolic  distribution  of  pressure 

A  force  F  applied  normally  to  the  contact  area  may  be  distributed 

as  pressure  p  given  by: 

2  2  4 

p  =  2F(a  -  p  )/ira  for  p  <  a 


0  for  p  >  a 


(3) 


which  resembles  a  paraboloid  of  revolution,  with  value  2F/rra  on  the  z 
axis,  zero  at  the  contact  edge.  This  case  was  mentioned  by  Boussinesq 
(1885  pl49)  who  calculated  the  normal  displacement  w  of  the  compressed 


surface  as 
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where  his  (  A f  2p)/u(  A  +  p)  =  2(1  -  v)/G.  Further  details  of  this 

field  do  not  appear  to  have  been  published,  but  may  readily  be  found  by 
using  Boussinesq*  harmonic  functions  as  described  by  Love  (1929). 

If  a  function  X  is  defined  by  the  equation 


Xdp , z)  = 


f" 

J 


p  log  (z  +  r!  dS 


(5) 


where  dS  is  an  element  of  the  contact  area  0  <  p  <  a,  r  its  distance 
from  the  field  point  [p ,  (p ,  z)  ,  and  if 
V(p,z)  =  dX/bz 
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p  r  1  dS 


(6) 


then  and  V  are  harmonic  functions  from  which  we  may  calculate  the 
stresses  and  displacements,  throughout  the  region  z  >  0,  caused  by  the 
given  applied  pressure  p.  The  stresses  are  given  by: 
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The  other  shear  stresses  are  zero.  For  the  displacements  u  and  w  we 


4itG  w 


2(1  -  v)V  - 


and  from  these  equations  the  whole  field  may  in  principle  be  determined,  by 

substituting  a  particular  form  such  as  (3)  for  the  pressure  p  in  (5) . 

But  it  is  not  always  easy  to  calculate  y  and  V  from  (5)  and  (6)  as 

explicit  functions  of  p  and  z,  although  V  has  been  determined  for  the 

3V 

Hertz  pressure  distribution,  and  found  for  the  case  of  a  unirorm 
pressure  (love,  1929).  The  function  V  is  also  known  for  pressures  due  to 
smooth  rigid  indenters  of  flat  or  conical  form  (Love  1939,  Green  1947)  but 
attempts  to  solve  the  present  case  of  parabolic  pressure  (3)  have  not  yet 
been  successful. 

However,  much  can  be  done  with  the  surface  and  axial  values  of  the 
functions,  which  are  readily  determined.  We  have,  on  p  =  0, 


V (o, z)  =  2irj  E“^ — ~h  Pi  dP 
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since  V  and  X  have  axial  symmetry  ani  satisfy  Laplace's  equation.  For 
the  function  }(  we  may  integrate  V  with  respect  to  z: 

?C(o,z)  =  V  dz  (10) 

F  F  2  2  3 /t>  2  2  2  t  4 

=  —  2z  (a  +  z  )  '  +  3a  z (a^  *  z  )'  -  2z 


4  2  2  t  2  2 

+  3a  log  (z  +  (a  +  z  )  *)  -  6  a  z  +  const. 


Equations  (7)  to  (10)  determine  the  stress  components  and  w  on  the  z  axis. 

To  find  their  values  at  other  points  we  observe  that  V  as  given  by 
(9) ,  but  regarded  as  a  function  of  a  complex  variable  z,  is  an  analytic 
function  of  z  in  a  domain  which  includes  the  origin.  It  is  therefore 
possible  to  express  V  as  an  integral  of  the  same  function  of 
(z  +  ip  cos  9)  (see  Whittaker  and  Watson,  1946,  p399) .  For,  if 
V ( 0 , z )  =  g(z)  in  (9) 


V(p,z)  =  —  |  g (z  +  ip  cos  9)  dc  (11) 

t. 

J  o 

in  a  region  surrounding  the  origin.  Similar  relations  may.  be  written 
3  V  3^V 

down  for  —  and  — since  they  also  are  harmonic,  but  this  does  not 

dZ  dZ”* 

apply  to  the  derivatives  with  respect  to  p. 

Although  these  integrals  may  be  cumbersome  for  an  arbitrary  point 


(p,z),  they  can  readily  be  evaluated  on  the  surface,  where  z  =  0  and  p  <  a 
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(in  agreement  with  Boussinesq  p.150) 
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where  E,  K  are  complete  elliptic  ir.t  ^rals  of  modulus  k  -  ft/ a.  Similarly 
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These  equations,  with  (7)  and  (8),  give  the  displacements  and 
stresses  on  the  contact  area.  At  points  [^q,o )  on  the  free  surface  outside 
the  contact  region  the  stress  components  and  Uj  take  their  usual  values 
corresponding  to  the  load  F  (see  Way,  1940) ,  and  w  is  found  from 
Boussinesq1  integral  (1885  p.1505. 

The  results  for  the  parabolic  pressure  distribution  are  summarised 

below. 

(i)  On  the  z  axis,  where  p  =  o: 
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as  for  any  symmetric  distribution  of  the  load  F  (Way  1940)  but  w  is  here: 
w  =  3(1  -  v)Fo  \  (2  -  3k2)  (1  -  k2)K  2(2k2  -  1)e1 /97i2Ga2k2 

where  the  elliptic  integrals  now  have  modulus  k  =  a/p. 

From  equation  (14^)  the  curvature  of  the  displaced  surface  at  the 

axis  is 

1/R  =  2(1  -  v)F/irGa3 

P 

and  the  vertical  displacement  there  is 

w  =  4(1  -  v)  F/3-irGa  (15) 

o 

so  that  this  field  ^gives  a  slightly  deeper  and  more  pointed  indent  them 
the  Hertz  formulae.  The  maximum  pressure,  at  the  centre  of  the  contact 
area,  is 

-  oz  =  2F  /tt  a 2  (16) 

or  twice  the  mean  pressure. 

This  completes  our  description  of  the  parabolic  pressure  field. 


which  is  used  as  a  correction  to  the  Hertz  solution  in  55.  But  first 


since  other  contact  problems  may  require  similar  corrections,  we  briefly 
describe  in  §4  a  few  simple  pressure  fields  and  some  useful  relations 
between  them.  This  section  may,  however,  be  omitted  by  the  reader  who 
wishes  to  proceed  to  §5. 

4.  Some  simple  pressure  fields 

In  linear  elasticity  theory  many  problems  are  solved  by  combining 
separate  known  stress  fields  in  suitable  proportions.  This  present  paper 
makes  use  of  the  Hertz  solution  and  the  parabolic  pressure  field,  but  there 
are  three  other  well-known  axially-symmetric  cases,  namely  the  flat  punch 
(Boussinesq  1885  p.158),  the  conical  punch  (Love  1939)  and  the  case  of 
uniform  pressure  (Boussinesq  p.t40,  Love  1929)  .  To  these  we  may  add  the 
linear  or  conical  pressure  field.  Some  interesting  properties  of  these 
fields  and  some  relations  between  them  are  worth  mentioning  here. 

Arranged  in  two  groups  of  three,  the  six  cases  are  displayed  in 
Table  I,  with  diagrams  of  the  pressure  distribution  and  the  shape  of  each 
indentation,  and  with  the  corresponding  formulae  alongside,  Since  these 
fields  may  be  combined  linearly  it  is  often  possible  to  build  the  shape  of 
a  required  pressure  distribution  or  given  indentation  from  these  known  fields 
The  process  is  simple  if  the  contact  area  is  the  same  in  each  case. 

Alternatively,  a  given  pressure  distribution  may  be  represented 
approximately  by  applying  uniform  pressure  (B^  in  Table  I)  to  smaller 
contact  areas  of  radii  a/n,  2a/n  ...  (n-l)a/n,  a,  so  making  a  stepwise 
approach  to  the  required  curve,  becoming  closer  as  n  increases.  Both 
and  (in  Table  Z)  may  be  regarded  as  integrals  of  B^  in  this  way,  but  it 
is  not  so  easy  to  build  a  required  indentation  shape,  such  as  a  parabola, 

A^,  or  cone,  A^,  by  superposing  solutions  for  varying  contact,  because 


the  non-zero  displacements  w  outside  the  contact  must  be  included  in  the 


10. 


calculation. 

Now  having  seen  the  parabolic  pressure  as  an  integral  of  simple 
uniform  ones  with  respect  to  a,  we  may  reverse  this  process  and  regard  the 
whole  stress  field  as  a  derivative  of  B^  or  B ^ ,  and  the  flat  punch  Aj  as 


a  derivative  of  Hertz,  as  follows. 

As  described  by  Love  (1929)  the  elastic  stress  field  of  a  particular 
pressure  distribution  may  be  derived  from  a  single  harmonic  function  x*  by 
using  eqns.  (5)  to  (8)  in  53.  For  the  Hertz  case. 


XH(p,z)  =  1  3  \2  (a2  -  Pj2)12  log  (z  +  r)  Pj  di  dPj 

2ita  I 

:/  o Jo 


where  r  is  the  distance  between  the  field  point  (p,  o,  z)  and  a  point 


(Pj»  <fr#  o)  in  the  contact  area.  If  we  write  this: 


3  3F  f  3  f 2V  .  2  2,  h  .  ,  .  . .  . 

3  XH  =  27  o  o  U  "  P1  )  log  lz  +  r)pi  d*  * 


1 


1 


then  a3xfl  is  a  continuous  function  of  a  for  any  j-fxec/  point  (p,z).  We  may 
therefore  differentiate  this  function  with  respect  to  a,  and  so  obtain  another 
harmonic  function  of  p,zj 


a 

3a 


,  3  .  3Fa  f 

(a  V  =  17  J 


2ir 


lo9<z  +  r)  .  . 

2  Th  P1  d*  d  P1 

o  Jo  (a2  -  Pj  V  1  1 


(since  the  pressure  is  zero  at  p^  =  a).  But  for  a  load  F  applied  with  a 
rigid  flat  circular  punch  (case  Aj)  the  fundamental  harmonic  function  is  x«» 


where 


f  a  f  2ir  log  _(Z  +  r) 

I  1  ~2  T7  P1  d*  dpi 

J  O  J  o  (a  -  Pj  ) 


F 

Pf2’ 

2ira 

«• 

9  (a3 

xH) 

3a 


(17) 


9a 


This  is  the  simple  relation  between  the  two  elastic  fields,  and  by 
interchanging  the  order  of  differentiation  it  is  possible  to  regard  a 
component  of  the  flat  punch  field  as  a  simple  derivative  of  the  corresponding 
Hertz  component. 

For  example,  the  displacement  in  the  contact  area,  for  the  flat 
punch  is  given  by: 


u  (punch)  =  — ~r  —  (a3  u  (Hertz) ) 

p  3a2  3a  p 


(1  -  2v)F  _3_  3 

,  2  .  _  3a  3 

3a  4irGp 


,  2  2  3/2. 

(a  -  p  )  ) 


(1  -  2v)F  ,  ,2  2S 

la  -  (a  ‘  p  >  1 


Similarly  the  function  V,  that  is  ,  for  the  flat  punch  is 


V  *  -  tan  1  (— ) 

Fa  z 


on  the  z  axis,  where  the  Hertz  function  is 


3F  2  2  -1 

V  =  — -  (  (a  +  z  )  tan  (a/2)  -  az) 

H  2a3 


So  again  V 


3a 


1  3  ,  3  „  , 

2  3a  (a  V 


These  two  related  fields  have  simple  algebraic  expressions  defining 
both  the  vertical  displacement  and  applied  pressure  in  the  contact  area, 
and  were  described  by  Boussinesq  as  the  two  simplest  cases  (1885  footnote 
p  206).  He  did  not  mention  the  case  of  the  conical  punch  (Love  (1939) 
Harding  and  Sneddon  (1945)  Green  (1949))  which  has  a  simple  linear 


displacement 


12. 


(1  -  v)  F 
_  2 


,  173  \ 

(  —  -  p) 


in  the  contact  circle,  perhaps  because  of  an  infinity  in  the  pressure. 


p  =  ~^~2  (log  (a  +  (a2  -  p'  )S  -  log  p) 

ira 


This  conical  punch  field  (A^  on  Table  I)  is  also  related  to  the  flat  punch 
solution  by  a  simple  differentiation 

XF  ‘  27  Ti  ^  *c>  1181 


where  xc  la  the  fundamental  function  for  the  conical  punch.  Equations 
similar  to  (18)  relate  corresponding  components  of  the  two  fields. 

The  three  solutions  Aj,  A^  and  of  Table  I  therefore  form  one  group, 
and  the  relations  (17)  and  (18)  between  them  can  be  used  to  reduce  tedious 
calculations,  or  to  check  a  derived  formula. 

The  other  three  solutions  and  B^  are  related  in  a  similar  way. 

For  if  Xp  is  the  fundamental  function  for  the  parabolic  pressure  problem 


(B2) ,  then 


3,4, 

T  Sa  <a  y 


is  the  function  for  Love's  case  of  uniform  pressure  (B^),  and  moreover 

■  rj  £  V  (,9) 


where  xT  relates  to  the  linear  (or  conical)  pressure  distribution  (B.) . 

The  related  stress  fields  B  j ,  B2  and  B^  of  this  second  family  do  not 
have  as  simple  expressions  for  the  displacement  as  the  A  group.  Instead, 
w  is  expressed  in  terms  of  complete  elliptic  integrals  of  modulus  (p/a) , 
and  the  calculation  of  stresses  is  more  difficult  than  in  the  cases  A. 


However,  the  axial  values  are  guite  tractable,  as  was  shown  In  S3 


be  found  from  V(0,z)  using  equations  (7).  For  an  arbitrary  point  (p,z) 
the  fields  of  3^  and  require  furtler  simplification  of  the  integrals. 

The  relations  (17),  (18)  and  (19)  between  these  different  fields  of 
elastic  contact  can  be  useful  in  practical  calculations  both  as  a  check 
and  as  a  wider  view.  The  parabolic  pressure  field  which  we  are  using  here 
can  now  be  seen  not  as  a  strange  new  field  but  in  its  rightful  place, 
being  related  to  the  uniform  pressure  problem  jSSBfc  as  Hertz  is  related  to 
the  flat  punch. 


5.  The  combined  solution 

This  section  may  be  read  as  if  directly  following  §3. 

A  linear  combination  of  the  two  elastic  fields  is  formed: 

(1  +  d)  Hertz  -  d  (parabolic)  (20) 

with  d  a  positive  parameter  which  is  to  be  determined.  Each  displacement, 
strain  or  stress  component  of  the  combined  field  is  of  the  form  (20) ,  with 
the  same  load  F  and  contact  radius  a  for  the  combined  as  for  the  separate 
fields. 

Then  if  the  combined  solution  is  to  give  an  indent  which  fits  against 
a  rigid  sphere  in  0  <  p  <  a,  the  best  value  of  d  must  be  found  for  each  F, 
that  is,  the  d  value  which  makes  the  displaced  surface  closest  to  a  sphere 
The  criterion  chosen  for  this  is  quite  simple.  Taking  x,  y  axes  not 
from  0  but  frcm  the  lowest  point  of  the  displaced  surface  in  fig.l,  the 
coordinates  of  a  surface  element  originally  at  (p,0)  are: 
x  «  p  ♦  u 

p  (21) 

y  -  »0  -  * 

where  w  is  the  vertical  displacement  at  p  -  0.  The  set  of  points  (x,y) 


corresponding  to  values  of  p  less  than  a  all  lie  on  a  sphere  of  radius  R 
with  centre  on  the  z  axis  if: 


2  2  2 
x^  +  (R  -  y)  =  R 

This  requires 

2  2 

(x  +  y  )/2y  =  constant  =  R  (22) 

and  this  has  been  used  to  find  d  for  various  values  of  the  load.  The 
fitting  process  determines  not  only  d  for  that  load  but  also  the  radius  R 
of  the  sphere  which  would  give  that  a  at  the  given  load  F.  For  very  small 
loads  d  is  zero,  since  the  Hertz  solution  is  accurate  there.  For  larger 
loads  d  increases  almost  linearly  with  the  mean  contact  pressure 4/ 
Values  of  d  have  been  found  in  this  way  at  various  loads  for  three 
values  of  Poisson's  ratio,  v  =  0.1.  0.25  and  0.4.  The  results  are  shown  in 
Table  II,  and  also  the  values  of  the  TCAlioa/'R  of  the  indent,  with  the  Hertz 
Ts&p  cl/Rh  for  comparison.  R  differs  from  R^  for  many  of  the  wider  contacts 
but  does  not  vary  in  a  simple  predictable  way.  A  numerical  correction 
factor  D(d)  is  defined  by  R  =  Rj^/Dfd),  and  this  also  is  shown  in  Table  II  # 
f  i  y  4  . 

The  same  results  are  displayed  graphically  in  fig. 3,  where  the 
straight  dashed  line  indicates  the  simple  proportionality  of  mean  pressure 
to  a/RH  as  in  eqn.  (2) ,  and  curves  are  drawn  for  the  modified  relations 
a/R  =  D (d) .  3(1  -  v)F/8  Ga2  (23) 

-  D(d) f/RH 

showing  the  effect  of  the  correction  factor  D(d).  This  equation  may  also 
be  written: 

a3  =  0(d)  3(1  -  v)F  R/8G 

-  0(d)  a3H  (24) 

showing  that  (D(d))1^3  is  the  correction  factor  for  a  as  calculated  from 


$ 

) 


known  F  and  R  by  the  Hertz  formula.  The  curves  are  drawn  up  to  a/R  *  0.8, 
although  not  many  materials  are  car  ble  of  perfect  elasticity  to  this 
extent. 

It  may  be  seen  from  fig. 3  that  for  v  =  0.1  the  ratio  a/R  is  always 
greater  than  the  Hertz  value,  while  for  v  =  0.4  it  is  less.  For  v  =  .25 
the  difference  changes  sign  at  about  a/R  =  0.5,  with  a/R  less  than  the 
Hertz  value  in  the  upper  range.  The  difference  in  the  lower  range  is  not 
distinguishable  in  the  figure,  but  a/R  exceeds  the  Hertz  value  by  about  2% 
at  a/R  'v  0.3. 

Although  D(d)  =  1  at  a/R  ^  .5,  so  that  a  =  a„  and  R  =  R„,  yet  there 
are  corrections  to  be  made  nevertheless.  For  at  this  point  d  llAS  T At 

value  0.27,  which  modifies  each  component  of  the  field  according  to  the 
formula  (20) .  For  example  the  maximum  pressure  under  the  indenter  is  now 
given  by 


(1  +  d) 


which  is  less  than  the  Hertz  value  by  9%,  and  the  depth  of  the  indentation, 
from  (1)  and  (f^*) ,  is  reduced  by  nearly  4%.  . 

The  values  of  d  in  Table  ZZ  and  fig. 4  show  that  the  stress  field  of 
Hertz  is  increasingly  modified  as  the  contact  area  widens,  requiring 
larger  proportions  of  the  parabolic  field  to  give  a  good  fit  to  the  sphere. 
But  while  d  is  increasing,  the  correcting  factor  D(d)  may  vary  from  unity 
in  either  direction.  This  means  that,  for  a  given  series  of  increasing 
loads,  a/R  becomes  greater  than  predicted  by  Hertz  if  v  is  small,  but  less 
than  the  Hertz  value  for  large  v,  and  changes  over  at  a/R  *v>  .5  for  v  *  .25. 

The  corrected  values  are  not  perfect  because  we  are  combining  only 
two  separate  curves  to  fit  a  circle.  After  finding  the  best  value  of  d, 
and  the  corresponding  R,  the  remaining  gap  or  overlap  of  sphere  and  surface 


is  given  by 


[x2  +  (R  -  y} 2]  -  R 

This  remaining  misfit  is  nearly  .0003R  at  a/R  =  .5,  and  .OOlR  at 
a/R  =  0.8.  Comparison  with  fig. 2  shows  the  degree  of  improvement. 

The  results  presented  in  Table  II  and  figs.  3,4  may  now  be  used  to 
estimate  errors  in  practical  cases.  An  example  is  demonstrated  in  the 
next  section. 

6.  An  example 

Let  us  suppose  that  the  elastic  modulus  E  of  a  pliant  material  is  to 
be  determined  experimentally  using  a  hard  spherical  indenter . 

From  the  Hertz  equations  (1)  and  (2)  we  have,  since  G  =  E/2(l  +  v) , 
wq  =  3(1  -  v2)F/4Ea 

=  a2/R  (25) 

and  so 

E  =  3(1  -  v2)F/4r‘!  w  3/2  (26) 

H  o 

where  wq  is  the  displacement  of  the  rigid  indenter  from  a  position  just 

touching  the  original  flat  surface  of  the  specimen  (see  Cousins,  Armstrong 

and  Robinson  1975) .  A  series  of  accurate  measurements  of  w  and  the 

o 

corresponding  load  F,  for  a  known  indenter  radius  R,  gives  us  values  of 
2 

E„/(l  •  v  ) ,  or  E  itself  if  v  is  known. 

H  p 

If  the  values  of  E^  so  determined  appear  to  rise  (or  fall)  with 
increasing  load,  then  it  is  of  interest  to  decide  whether  this  is  due  to 
anelastic  behaviour  of  the  material  or  merely  to  the  inaccuracy  of  the 
Hertz  formulae.  This  question  is  readily  answered  using  fig. 3,  which 
displays  the  calculated  curves  (23),  and  the  corresponding  values  of  d  rVi  Utj  4. 
Consider  first  the  curve  for  v  -  0.1.  The  factor  D(d)  is  unity  for 
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small  contact,  increases  to  1.05  at  a/R  ^  .5  where  d  =  .3,  then  falls 
again  to  1.01  at  a/R  ^  .8,  d  =  .55.  For  w  we  have  the  combined  form  (20), 

w  =  (1  +  d)  w  (Hertz)  -  d  w  (parabolic) 

o  o  o 

=  (1  -  0. 1317d)  3(1  -  v2)F/4Ea 

«  (1  -  0. 1317d)  a2/RD(d) 

Substituting  a  =  [D(d)  R  wq/(.1  -  0.1317d)]  ^ 

gives  an  equation  for  E  in  the  corrected  form: 

E  =  (1  -  0.1317d)3/2  3(1  -  m^F^R*5  wq3/2  Did)**) 

=  (1  -  0.1317d)3/2  E  /D (d) 15  (28) 

□ 

From  this  it  appears  that,  for  v  =  0.1,  the  "measured"  values  of  E  , 

n 

obtained  by  using  eqn  (26),  would  exceed  the  true  values  E  as  the  load 
increased,  the  error  being  9%  for  a/R  ^  .5  and  12%  for  a/R  ^  .8. 

Repeating  these  calculations  for  v  =  .25  we  find  similar  behaviour  but 
a  smaller  effect,  errors  now  6%  and  8%.  Finally  for  v  =  .4  the  error 
reaches  only  1.4%  at  a/R  *\<  .5  then  decreases  to  1%  at  a/R  ^  .8.  The  Hertz 
formula  therefore  seems  reliable  for  E  determination  when "v  is  large,  and 
this  is  fortunate,  since  such  v  occur  in  some  very  pliant  glassy  polymers 
for  which  an  elastic  indentation  may  well  reach  high  values  of  a/R.  The 
error  remains  small  because  the  corrections  in  numerator  and  denominator  of 
eqn  (28)  are  of  similar  magnitude  in  this  case. 
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7.  Conclusion 

The  classic  formulae  of  Hertz,  applied  to  the  indentation  of  an 
elastic  half-space  by  a  rigid  sphere,  have  been  modified  to  give  greater 
accuracy  for  contact  areas  beyond  the  usual  range.  For  this  purpose  a 
large  part  of  the  elastic  field  of  a  parabolic  pressure  distribution  has 
been  described  in  detail,  and  also  discussed  more  generally  in  relation  to 
other  contact  problems. 

The  parabolic  field  has  been  subtracted  in  small  proportions  from  the 
Hertz  field,  making  use  of  the  linearity  of  the  system.  The  ratios 
required  to  give  accurate  fit  to  the  sphere  at  wide  contacts  have  been 
calculated  for  three  different  values  of  Poisson's  ratio,  and  displayed 
graphically  for  practical  application.  Comparing  the  modified  solution 
with  that  of  Hertz,  the  changes  in  the  radius  of  contact  are  in  general 
small,  with  maximum  4%,  but  larger  corrections  were  found  for  the  maximum 
pressure  and  for  the  indentation  depth. 

As  modern  pliant  materials  extend  the  elastic  range  of  indentation, 
and  modern  methods  increase  the  accuracy  of  measurement,  these  results  may 
be  useful  for  comparison  or  reassurance.  The  theory  is  first  order  linear 
elasticity  throughout. 
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Fig.  2 


Fig.  3 


Fig.  4 
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Diagram  of  indented  surface  showing  axes  and  extent  of  contact. 

The  errors  in  the  Hertz  theory  as  a/R  increases.  Positive 
values  of  R  -  R  indicate  overlap,  negative  ones  a  gap. 

Mean  pressure  vs  modified  extent  of  contact  a/R. 

Variation  of  d  with  mean  pressure  for  various  v. 


*0 


•M  V 


>  *V'  ' 


IA  • 


■Wi 


/- 

Sl'Vi,  ’ 

fe  tut  ■ 


.-(>■ 


m  ; 


V“4«  '  v’  ,  • 

,'  r  *;  *•  '  •;:>->•  •  .  ■■• 

•  :M.  ’•/• 

^  -  t\  y, .  ,  V  ’  ■■*  .  i  *’  * 

'  *  :  A.’  ’-J’-  " 

1  'V.-V’?;.  «'.»»•.?  *  S 


J'.f5  ' 


%  ■ 


m  •  "•;  >.'  • 1  • . *,...  . 

,w  ‘  ■  ;•■•  «  ■*■*■$*«*  -v  * .  « ■•  ;.  .v*v ,  : 

*!;,  -**  ^  ■  \  a-  %  %  .-/■  .  a  ■  a’-'.. 

.  A  •  -  ■*  ■  .1  _  **« , v  '  J  > 1  ..  t 


■i 


4  t'-i!’  /(..•  ;  >  *r  *  ;  "A 

,  1  a 


.;•,•*  'Vu  ‘Ji  •  •*’&' Jr'  'lj<  A  tfr/  Am-:' 


.  .?$$**•. ^ >*3;  ■  i5:;: 

*,  :•  .»»,  •/>  ^  <»h;>  ,v  ^  ^ .  ■ 


*  ,vfr-»  •■  ■  *  ■vm*  ■  .-1 


» 

t, *  K 

:*r7  /t  •■'  ;  ,  > 
,  «  ■  !<., 

■•*.  7>  ^  ^  r*  _ 


>.i  .»•'.*  • 


